Magnetic and transport properties of a one dimensional frustrated t-J model for 

vanadate nanotubes 
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We propose a one-dimensional model consisting of a chain with a t-J Hamiltonian coupled to a 
Heisenberg chain in a frustrated geometry to describe the appearance of the ferromagnetic phase 
which has been experimentally observed in vanadate nanotubes. This model contains a mechanism 
of frustration suppressed by doping suggested by L. Krusin-Elbaum, et al. [Nature 431, 672 (2004)]. 
We study, using numerical techniques in small clusters, the relation between magnetic order and 
transport properties in the proposed model, and we perform a detailed comparison of the physical 
properties of this model with those of the ferromagnetic Kondo lattice model. For this comparison, 
a number of results for the latter model, obtained using the same numerical techniques, will be 
provided to complement those results already available in the literature. We conclude that it does 
not appear to be a true ferromagnetic order in the proposed model, but rather an incommensu- 
rate ferrimagnetic one, and contrary to what happens in the ferromagnetic Kondo lattice model, 
electronic transport is somewhat suppressed by this ferrimagnetic order. 

PACS numbers: 71.10.-w, 71.27.-fa, 75.30.Kz, 75.40.Mg 



I. INTRODUCTION 

The possibility of taking advantage of the spin of the 
electrons in addition to their charge in order to develop 
devices vi^ith new capabilities, which is exploited in the 
emergent field of spintronics^, has in turn given impulse 
to the search for new materials with metallic ferromag- 
netic (FM) phases. 

An extensively studied family of compounds which 
presents these properties is the family of the so called 
"manganites" , for example Lai_a;Caa:Mn03i^ This fam- 
ily of manganese oxides has recently received additional 
interest for its property of "colossal" magnetoresistance.^ 

Taking into account the vast diversity of physical be- 
havior found in transition metal oxides, it seems natural 
to look for this kind of properties within these materi- 
als. Following this program, it was found that a vana- 
dium oxide {yOx) presents upon doping a ferromagnetic 
phase at room temperature^ This compound, as well as 
manganiteS}^ has been synthesized as very thin cylinders 
or nanowires thus opening the way to their use as ferro- 
magnetic leads in nanoscopic spin valves^ among other 
devices. In this compound, double-layered nanotube 
"walls" of vanadium oxide are separated by dodecylam- 
monium chains i^i^ii^ This structure is reminiscent of the 
carbon nanotubes which are the essential components in 
important devices that have already been developed. 

In manganites, it is now widely accepted that most of 
their relevant physical properties, particularly the prop- 
erty of magnetoresistance, can be explained by a single- 
orbital ferromagnetic Kondo lattice (FKL) model, ^^'^"^ 
which can in turn lead to the well-known "double ex- 
change" (DE) model in the limit of infinitely large Hund 
coupling. Of course, it is necessary to resort to at least 
a two-orbital model to describe properties related to or- 
bital orderingjii 



The purpose of the present work is to examine an al- 
ternative mechanism to the one contained in the FKL or 
DE models. This mechanism was suggested in Ref. and 
is essentially based in the suppression of frustrated anti- 
ferromagnetic couplings in VO^; walls by electron or hole 
doping. In order to study this mechanism we propose a 
model inspired by the VO^: compound but we would like 
to emphasize that this model is not aimed at a faithful 
description of this material. 

Our model, shown in Fig. [1] below, is a simplified one- 
dimensional (ID) version of the VO^; structure. ^'^'^^ This 
model differs from a realistic model for VO^ not only in 
its dimensionality but also in a number of details such 
as the zig-zag geometry of the chains, the alternation be- 
tween Vi and V2 along the chain, and a possible dimer- 
ization which may explain the spin gap experimentally 
observed in this material. Nevertheless, we believe this 
model will capture the essential features of the mecha- 
nism suggested in Ref. and at least qualitatively some 
properties of the real materials. In fact, V1-V2 chains 
are structures present in VO^: although coupled between 
them in a complex way. In support of the present ap- 
proach, we would like to recall that results for the one- 
dimensional FKL modeiJ'^ have reproduced many essen- 
tial properties observed in manganites. 

In the present manuscript we will mainly address, us- 
ing numerical techniques in small clusters, the relation 
between magnetic order and transport properties in the 
proposed model, and we will perform a detailed compar- 
ison of the physical properties of the present model with 
those of the FKL model. For this comparison, a number 
of results for the latter model, obtained with the same 
numerical techniques, will be provided to complement 
those results already available in the literature.— i^^ii^ 

The paper is organized as follows. In Section we 
present the model studied and we provide details of the 



2 



v V 

FIG. 1: A ID model for VO^. The upper chain, with Vi 
and V2 ions, is described by a t — J model, while the lower 
chain, with V2 and V3 ions, is described by an AF Heisenberg 
model. 

numerical techniques employed. In Section IIIII we show 
the existence of various magnetic phases in the model 
studied. Transport properties, in particular optical con- 
ductivity, are examined in Section HVl Section |V] is dedi- 
cated to emphasize differences between the present model 
and the FKL model. 



II. MODEL AND METHODS 

The ingredients of our effective one-band model for 
vanadium orbitals only, which are common to VO^;, are 
the following (see Fig. [T|). In the first place, the sites of 
the system correspond to three types of vanadium ions, 
Vi, V2, and V3, with a hierarchy of hole on-site energies, 
ei < £2 < es."* Vi ions are then the most likely to be 
occupied by holes upon doping. 

The couplings between vanadium ions are mediated 
by oxygens and the relative values are determined by 
the quasi-pyramidal (Vi, V2) and tetrahedral (V3) co- 
ordination between vanadium and oxygens. In the real 
YOx, Vi and V2 ions form zig-zag chains running along 
a given direction in one layer and in the perpendicular 
direction in the nearest neighbor layers. These layers of 
chains are connected through V3 ions which are located 
between them. The couplings between Vi, V2 and V3 
form a frustrated triangle of antiferromagnetic (AF) in- 
teractions. 

In order to simplify the calculations we assume that 
double-occupancy is forbidden, which is reasonable for 
d-orbitals in transition metals, particularly close to half- 
filling. Then, the model for V-ions, with the geometry of 
Fig. [l]is defined by the Hamiltonian: 

i.<7 i 

s,2/-l • Ss^2i + Ss^2/ ■ Ss,2l + l) 

I 

+ E('/lSs,2/-l ■ 82;-! + J2^s,2l-1 ■ S21) 
I 

+ ^e,(l-n,) (1) 

i 

where cj^. and Cia are electron creation and annihilation 



operators with the constraint of no double occupancy on 
Vi and V2 sites, and Sj (Ssj) are the spin-1/2 opera- 
tors on the t-J (Heisenberg) chain. The first two terms 
correspond to the t — J Hamiltonian for the upper chain 
in Fig. [1] the third term to the Heisenberg Hamiltonian 
for the lower chain with odd (even) numbered sites corre- 
sponding to V3 (V2) sites, and the forth term to the ex- 
change couplings between these two chains. For simplic- 
ity, in the present study we will not consider dimerization 
in the model, and hence we adopt the hopping integrals 
ti = <2 = i and the exchange couplings Ki = K2 = K , 
and Ji = J2. Besides, from the structure of VO^;, it 
should be J = J2 as well. All exchange couplings, K 
and J, are assumed antiferromagnetic. K — \ \s adopted 
as the scale of energies. e2 = gives the reference for 
the on-site energies; ei < 0, which will be a variable pa- 
rameter in our study, and — 00, implying a forbidden 
hole occupancy on V3 sites. The £3 = 00 condition would 
cause localization of holes eventually occupying V2 sites 
on the lower chain in Fig. [T] and then it would lead to an 
effective cut of this chain. This is the reason why hole 
occupancy on V2 sites on the lower chain is forbidden in 
our model. In VOx, these V2 sites would correspond to 
V1-V2 chains running in a perpendicular direction of the 
plane of Fig. [1] 

At half-filling, i.e. when all sites are single occupied, 
the system is reduced to a frustrated spin system. In the 
proposed model, which could be termed a frustrated t-J 
(FTJ) model, charge carriers are introduced by doping 
with holes the t — J chain. Most of the results reported 
below correspond to the doping fraction x = 0.4, where 
X is defined as the number of doped holes divided by the 
number of sites of each chain, L. Although most holes 
would go to Vi sites in YOx , in our model there would be 
a finite probability of the holes going on the V2 sites and 
the hopping V1-V2 would then be possible. In any case, 
the situation of localized holes on the Vi sites could be re- 
covered in the t = limit. The inclusion of the hopping 
term in our model makes it more general and eventu- 
ally its properties could be relevant for other compounds 
sharing similar structural features with VO^;. In our zero 
temperature study, the inclusion of this hopping is also a 
mimic of the activated transport that takes place in the 
real material. 

Model ([T]) was studied in 2 x L clusters by exact diago- 
nalization (ED), using the Lanczos algorithm, with peri- 
odic boundary conditions (PBC) and by density matrix 
renormalization group (DMRG)^ with open boundary 
conditions (OBC). Using ED we were able to compute 
static and dynamical properties on the 2 x 10 cluster, 
and using DMRG we studied 2 x 20 and 2 x 40 clusters 
although in this case only static properties were consid- 
ered. In the results reported below obtained using the 
DMRG method we have retained 300-400 states in the 
truncation procedure. 

To compute the total spin S we have adopted two 
procedures, (i) compute directly S • S = S{S + l)/2, at 
= and (ii) compute S as the maximum value of 
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at which the ground state energy is recovered, which im- 
pUes diagonaUzing the Hamihonian in several subspaces 
of fixed S\ 

Spin-spin correlations {SjSq) are computed exactly 
with the Lanczos algorithm in the ground state, in 
the 5^ = subspace. The static magnetic struc- 
ture factor x(q) is the Fourier transform x(q) — 
V^E„(^|^f)exp(zq-(rj-ri)), (N = 2L) where rj 
(q) are the real space (reciprocal) vectors of a rectangular 
ladder of L x 2 sites where the lattice of Fig. [1] assuming 
equidistant sites, can be embedded. We have verified that 
the results reported below do not depend qualitatively of 
geometrical details (relative spatial position of the ions) 
of Fig. [T] Within DMRG, we measured spin-spin correla- 
tions in the = subspace from one of the two central 
sites of the cluster, one belonging to the t-J chain and the 
other to the spin chain. In this case the magnetic struc- 
ture factor is computed as x(q) = J2j i'^j'^o) ^^P ' J); 
where "0" is one of the central sites on the spin chain 
and j belongs to the spin chain. In order to compare 
with results for the FKL model, in most of the results 
reported below for this model we have restricted the sum 
to sites on the spin chain, but they are qualitatively sim- 
ilar to the ones obtained with the full Fourier transform 
involving all sites. 

A measure of the magnetic order on each chain is pro- 
vided by the sum of the correlations: 

Sh ^{SliSlg) (2) 

i 

where the sum extends over the sites of the spin (V3- 
V2) chain and q is the spin of a V3 site on this chain. 
Sh — corresponds to an AF ordering and Sh > to a 
FM or a ferrimagnetic ordering of the spin chain. 

An important property related to transport is the real 
part of the optical conductivity, (j{uj) = DS{uj)+a'^'^^ {ui), 
where a^'^^{uj) is defined as the spectral function: 

a-^(c.) = - ^ |(^„|j|^o)P'5(a; - {E„ - Eo)) (3) 

71 

where \ipn) are eigenvectors of the Hamiltonian ([1]) with 
energies En, and 

J ^ it^{cl+ia^^<y - H.c.) 

i 

is the current operator (in units of e = 1) and the sum 
extends over the t — J chain in ([1]). a'^'^^ [uj) was computed 
using the continued fraction formalismii^ 

From the optical conductivity, the Drude weight can 
be computed as!^>i^ 

- = Uek -- I (4) 
n N TT J 

where Ej^ is the kinetic energy, i.e., the negative ground 
state average of the hopping term of the Hamiltonian 
(first term in H])). Ek is expressed in units of K, and D 



in units of e^K. As usual in the continued fraction cal- 
culation of spectral functions a Lorentzian broadening of 
the discrete peaks was adopted. In all the results shown 
below we took a Lorentzian width 6 = 0.1. 

As mentioned above, when Ji = J2 = 0, model ([1]) re- 
duces to uncoupled t—J (where K plays the role of J) and 
Heisenberg chains. We are also interested in comparing 
our model with the FKL model with spin- 1/2 localized 
spins defined as: 

(5) 

where Jh > is the Hund coupling, are the localized 
spins and s.^ are the spins of the conduction electrons. 
In this case the filling n is defined as the filling of the 
conduction chain. Notice that the t2g electrons in man- 
ganites are usually modeled with spin-3/2 operatorsJ^ 
The frustrated t-J model is perhaps more resemblant 
to an extended version of model ^ which includes an 
AF Heisenberg exchange between localized spinsJ^'SSi'^ 
Model ID) would also be closer to a version of the FKL 
model which includes a Hubbard repulsion U on the con- 
duction orbitalSf22j^ since the t — J Hamiltonian appears 
at an effective level when U ^ t?"^ However, we would 
like to emphasize that the proposed FTJ model was not 
obtained from the U ^ t limit of an underlying Hubbard- 
like model. In fact, most of the results reported below 
belong to the case oi t < K. 

All the results reported below correspond to the model 
([T]) except otherwise stated. 

III. MAGNETIC PROPERTIES 

We first examine the magnetic properties of this model, 
in particular the presence of a ferromagnetic phase, as 
detected experimentally in Ref. 0. The suggested mech- 
anism that causes this FM phase works as follows. At 
half-filling, due to the frustrated AF interactions, the 
system would be an antiferromagnet since it reduces to 
relatively isolated spin chains (V1-V2 chains in the real 
materials), or to a gapped spin liquid if these chains are 
dimerized. Upon electron or hole doping, Vi sites become 
non-magnetic, the frustration in the V1-V2-V3 triangles 
disappears, and the AF order in the V2-V3 subsystem of 
Fig. [1] implies an excess of the z-component of the total 
spin along one direction thus leading to ferromagnetism 
or more properly, to a ferrimagnetic order. 

This argument is expressed in the atomic limit. For 
non-zero values of t, there is a competition between the 
onsite energy, which tries to localize holes, and the kinetic 
term which tries to delocalize them. In addition, as in any 
t — J-like model there is a competition between kinetic 
and magnetic energies. 

Let us first examine the undoped system. Notice that 
by adding interactions between the V2 sites in the bot- 
tom chain of Fig.[Tl and for the case J = K, one recovers 
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FIG. 2: (Color online) Spin-spin correlations (SjSo) for the 
undoped system, (a) along the V1-V2 chain; (b) along the 
V3-V2 chain (V2 reference site). Results were obtained for 
the 2 X 14 cluster with PBC (open circles) and for the 2 x 20 
cluster with OBC (filled circles), J/K = 0.3. (c) Spin-spin 
correlations between a V3 site and sites on the V1-V2 chain, 
2 X 14 cluster with PBC, J/K = 0.3 (squares), J/K = 1.0 
(triangles). The normalization (So'S'o)=l was adopted. 

the frustrated AF Heisenberg ladder which has been ex- 
tensively studiedJS. The ground state of this model is 
a singlet with a spin gap. If J < if , then we have 
an asymmetric frustrated ladder where we can expect 
a similar behavior. A singlet ground state could also be 
expected for our frustrated model. Some handwaving ar- 
guments can be put forward to support this in the case 
J < K. In this case, an AF order on the upper chain of 
Fig. [1] would dominate each frustrated triangle and the 
effective interactions with couplings Ji and J2 would be 
small. Then, both chains are essentially uncoupled and 
also the lower chain would present the typical AF or- 
dering of Heisenberg chains. Our numerical calculations 
support this picture. In Fig. [2] we show that for both 
clusters with PBC and OBC, J/K = 0.3, the upper and 
lower chains have the spin-spin correlations (S^Sq) corre- 
sponding to an AF order (Fig. [H^a) and (b) respectively) . 
In Fig. [2lc) we show that (SjSq) between a V3 site and 
sites of the upper chain are indeed relatively small. In 
fact, we have also determined that this behavior extends 
up to J/K = 1.0, as shown in Fig.[2)Jc). Notice the good 
agreement between results obtained for PBC (with exact 
diagonalization) and for OBC (with DMRG) where the 
reference site is one of the two central sites. It is worth to 
note that the previous argument would also imply that 
the undoped system is spin gapless and our numerical 
calculations seem to agree with this possibility. 

In Fig. [3] the phase diagram of model ^ in the plane 
{t, J} is shown for the 2 x 10 cluster and 4 holes (doping 
fraction x = 0.4). For ei = -2 (Fig. [Sja)), there is 
an abrupt crossover from the singlet state {S — 0) to 
the ferrimagnetic state with the maximum S expected at 
this doping, S = 2. As t is increased the holes tend to be 
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FIG. 3: (Color online) Phase diagram of model ([T| on the 
2 x 10 cluster with 4 holes, (a) ei = -2, (b) ei = -1. The 
inset shows Sh (defined by Eq. ([2])) for t = 0.3 (circles), 0.5 
(squares) and 0.7 (triangles) as a function of J. 



less localized on Vi sites, and then the crossover takes 
place at larger values of J. When |ei| is reduced, hole 
localization becomes less favourable, and the crossover 
takes place at even larger values of J for a given t, as 
can be seen in Fig. [Sljb), which corresponds to ei = —1. 
More interesting is the fact that, for small values of t, i < 
0.3, there is an intermediate region with S = 1 between 
the 5 = and S = 2 regions. The behavior of Sh for 
several values of t as a function of J is shown in the 
inset of Fig. [31 This quantity presents a jump at the 
5 = 0/5* = 2 crossover, and since the nearest neighbor 
correlation (SgiS^Q) < 0, then there is a ferrimagnetic 
ordering of the spin (V2-V3) chain. That is, a snapshot 
of this system would show an alignment of the spins on 
the spin chain and at the same time an alignment of the 
spins of the t-J chain with opposite orientation to the 
previous one. 

As the cluster size is increased, the intervening regions 
with < S < Smax develop a "staircase", as it can be 
seen in Fig. ID In Fig. HJa), the phase diagram in the 
plane {t, J} for the 2 x 20 cluster, with the same dop- 
ing as before, x = 0.4, is shown. These results were 
obtained with DMRG and the clusters have OBC along 
both directions. The crossovers in the phase diagram of 
the 2 x 20 cluster were determined by comparing the 
energies obtained for different S"^. In this model, as 
well as in the FKL model discussed below, DMRG has 
a slow convergence and there is a tendency for low 5^ 
subspaces to fall into metastable states. Then, in or- 
der to determine the value of S it is important to go to 
large 5*^ and then to extrapolate the difference in energy, 
/\E = EiS"") - E{S'' = 0), to zero. This procedure is 
illustrated in the inset in Fig. Ufa). In this way we have 
also computed the total spin for some parameter sets, in- 
dicated in Fig. mja), in the 2 x 40 cluster. These values 
of S for the 2 x 20 and 2 x 40 clusters, together with 
the previous one for the 2 x 10 cluster, seem to suggest 
that the ferrimagnetic phase will survive in the thermo- 
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FIG. 4: (Color online) (a) Schematic phase diagram of model 
(dl on the 2 X 20 cluster with 8 holes, ei = -1. At t ^ 0.6, 
J = 0.9 and 0.7 (points indicated with stars), 5 = 8 and 7 
respectively on the 2 x 40 cluster with 16 holes, ei = — 1. 
The inset shows AE (defined in the text) vs. Sz for t — 0.2, 
J = 0.2 (filled symbols) and J = 0.4 (open symbols), 2 x 20 
cluster, (b) Spin-spin correlations on the 2 x 40 cluster with 
16 holes, ei = —1, t = 0.6 and values of J indicated on the 
plot. The normalization {SqSq)=1 was adopted. The curves 
have been shifted for the sake of clarity. 
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FIG. 5: (Color online) Static magnetic structure factor xi'i) 
as a function of J along cuts of constant t on the phase dia- 
grams of Fig. 131 2 X 10 cluster, 4 holes, (a) ei = -2, t = 0.8, 
(b) ei = —1, t — 0.2. The values of q are indicated on the 
plot, in units of 7r/5. The inset shows xil) for the ferro- 
magnetic Kondo lattice model as a function of Jh, L = 10, 
n = 0.6, qx ~ S (circles), 2 (triangles) and 1 (diamonds) in 
units of tt/5. 



dynamic limit. It is interesting to note that we have also 
observed this staircase in exact diagonalization studies 
of the FKL model on the L = 10 cluster (see below). In 
the suggested mechanism for ferrimagnetism, the mag- 
netic order in the spin chain would still be an AF one. 
It can be seen in Fig.HJb) that the spin-spin correlations 
{SjSq) for the 2 x 40 clusters, x — 0.4, are indeed essen- 
tially AF except for the presence of a small kink. The 
same behavior was found for the 2x10 cluster with PBC 
and the same doping in the S = 2 region. 

In order to study more systematically the behavior of 
the spin-spin correlations in the regions of different total 
spin let us turn to study the static magnetic structure 
factor, x(q)- Results for in the S*^ = subspace 

along a "cut" at t = 0.8 in Fig. EJa), and at t = 0.2 in 
Fig. [S^b) are depicted in Fig. [5Ua) and (b) respectively. 
These results were obtained by exact diagonalization and 
the first definition of discussed in Section [ll] was 

adopted. In both cases considered, in the S = region 
(small values of J/K), x(q) is maximum at = tt [x is 
along the chain direction), indicating an AF order, irre- 
spective of gj, = and tt. In Fig.[51[a), in the S = 2 region 
{J/K > 0.7), the peak of x(q) shifts away from = tt, 
indicating an incommensurate (IC) order. In Fig. [5jb), 
x(7r, 0) has a sudden drop upon entering in the 5 = 1 
region {J/K ~ 0.3) while x(37r/5,0) and x(47r/5,0) in- 
crease. A similar behavior occurs when entering in the 
S = 2 region {J/K > 0.6) but now the maximum of 
x(q) is located at q = (47r/5, 0) corresponding to the 
IC order. This behavior do not depend qualitatively on 
other choices of the set of reciprocal vectors entering in 
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FIG. 6; (Color online) Static magnetic structure factor x{l) 
computed with DMRG (a) in the 2 x 40 cluster, for several 
values of J indicated on the plot, ei — —1, and t — 0.3; (b) 
in the 2 x 20 cluster for ei = —1, t = 0.4, J = 1 at several 
hole doping Nh indicated on the plot; (c) for the FKL model 
in the L = 20 chain, x — 0.6, for several Jh indicated on the 
plot, (d) Total spin of the ground state of the FKL model as 
a function oi Jh, L — 2Q, x = 0.6. 



the Fourier transform. The inset shows x{q) for local- 
ized spins for the FKL model as a function of Jh for the 
same chain length, L — IQ, and filling n = 0.6. In this 
case, q is reduced from 37r/5 to its minimum possible 
nonzero value, 7r/5, which corresponds to FM order in 
the 5^ = subspace with PBC. It is clear the presence 
of two crossovers, at Ju 6.5 and at 19.5, where also 
the total spin of the system jumps from to 1 and then 
from 1 to 2, respectively. 



Similar behavior of the static magnetic structure factor 
is observed in the 2 x 20 and 2 x 40 clusters with OBC. 
Fig. [6ja) shows x{q) in the 2 x 40 as a function of qx for 
various values of J, at t = 0.3, obtained with DMRG. In 
this case only the spins of the spin chain have been in- 
cluded in the Fourier transform but results involving all 
spin-spin correlations are similar. It can be seen that, for 
a fixed t, xil) remains approximately unchanged during 
some intervals of J. Each of these intervals of J can be 
related to a given value of the total spin S, but there 
is not necessarily a one-to-one correspondance. As hole 
doping increases, the IC peak of xil) becomes sharper 
but it remains at q = {{L ± l)7r/L,0), as it can be seen 
in Fig. [6l^b). We will term these states with S* > as IC 
ferrimagnetic states. This behavior, which was also ob- 
served for the other clusters studied, is different from the 
one reported for the ID FKL modelj^^ In the FKL model, 
when the electron occupation is reduced from half-filling 
at a fixed value of the Hund coupling, the peak of the 
static magnetic structure factor, located at 2kF — mr, 
is reduced from its AF value = tt to g^, = in the 
FM region. In Fig. ^c) we show for comparison the evo- 
lution of xil) for the FKL model as Jh is increased at 
an equivalent cluster size and filling. In contrast with 
Fig. El^a), the peak of x{q) moves away from its IC po- 
sition as it approaches the crossover to the FM phase 
located at Jh ~ 18. Inside this phase, x{q) presents 
the typical form of a FM order in the — subspace. 
In Fig. [nijd) we show the corresponding evolution of the 
total spin of the ground state. Although DMRG calcu- 
lations have a very slow convergence for this model, a 
problem that is also present in Lanczos diagonalization, 
there are apparently regions in Jh with intermediate val- 
ues of S which may correspond to the various positions 
of the peak of xio) shown in Fig. [S{c). However, this 
staircase in the values of S is not as pronounced as the 
one for model ([1]). Moreover, these results could be af- 
fected by finite size effects. These differences between 
the two models illustrate the IC ferrimagnetic character 
of the ground state in the present model as opposed to 
the true ferromagnetic state present in the FKL model. 



IV. TRANSPORT PROPERTIES 

For the possible applications of VO^; in spintronics de- 
vices, it is essential to determine the correlation between 
its magnetic and transport properties. That is, one would 
expect that, as in manganites, transport is enhanced as 
the total spin of the ground state increases. 

An elementary indication of charge mobility is pro- 
vided by the kinetic energy, Ek- Fig.[7]Ja) and (b) show 
Ek obtained by ED for the 2 x 10 cluster with 4 holes 
in the S*^ = sector for ei = —2 and ei = —1 respec- 
tively. It can be seen in Fig. [7] that as J is increased for 
a fixed t there are step-like reductions of Ek at values 
of J which exactly coincide with the boundaries between 
regions with increasing values of S in the phase diagrams 
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FIG. 7: (Color online) Kinetic energy (in units of K) as a func- 
tion of J on the 2 x 10 cluster with 4 holes, t = 0.1, 0.2, . . . , 1.0 
from bottom to top (a) ei — —2, (b) ei = —1. 

shown in Fig. [31 That is, the increase of S is accom- 
panied by an increased hole localization. Although it is 
not perceptible in the scale of Fig. [71 in addition to these 
steps, Ek slowly decreases as J is increased due to the 
competing magnetic energy. 

For small values of the hopping constant t {t < 0.2), 
small increases of Ek as a function of J can also be 
observed in Fig. [Tl^b). This behavior can be understood 
from the fact that in the regime t << K, in order to 
gain magnetic energy, the spins on the t-J chain have 
the tendency to form hole-free islands thus leading to a 
phase separated state. The exchange couplings Ji and J2 
connecting the t — J chain with the Heisenberg chain are 
competing with K and then opposite to the formation 
of those islands. Then, a,s Ji — J2 — J increase there 
will be at some point a breaking of that phase separated 
state thus favoring hole delocalization and an increase in 
kinetic energy. 

A similar behavior of the kinetic energy has been ob- 
served in the 2 x 20 cluster with 8 holes, and for the 
2 X 40 cluster with 16 holes as shown in Fig. [H^a) and (b) 
respectively. In order to clearly appreciate the evolution 
of Ek for various values of t, the values of Ek have been 
shifted by multiples of t as indicated in the vertical axis. 
For a given value of t there is a general decrease of Ek as 
J increases, with some sudden drops which coincide with 
the successive increases of S shown in the phase diagram 
of Fig. [H It can be seen in the inset of Fig. [H that the 
kinetic energy per doped hole also decreases as hole dop- 
ing is increased. In the FKL model, Ek also decreases as 
Jh increases at a given density^^ due again to the com- 
petition with the magnetic energy (see further discussion 
below) but Ek per electron increases for a fixed Jh- 

Fig. [HI shows the regular part of the optical conductiv- 
ity a^'^^ioj) computed on the 2 x 10 cluster with 4 holes, 
(a) t = 0.4, (b) t = 0.9, ei — —1. Excluding the peak at 
Lu = which should be replaced by the Drude weight 
a^'^^{ijj) and cr(w) coincide, so we will use (j{uj) in the 
following. For both values of (t{uj) presents a broad 
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FIG. 8; (Color online) Kinetic energy Ek (in units of K) as 
a function of J, for several values of t indicated on the plot, 
ei — —1, (a) for the 2 x 20 cluster with 8 holes, (b) for the 
2 X 40 cluster with 16 holes. Ek has been shifted as indicated 
in the label of the y axis for clarity. The inset shows Ek per 
hole in the 2 x 20 cluster as a function of the number of holes. 
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FIG. 10: (Color online) (a) Drude weight (open symbols) and 
the integral of a^''^{u}) (full symbols) for model ([T]) for the 
2 x 10 cluster with 4 holes, ei = —1, as a function of J, t — 0.4 
(circles) and t = 0.9 (squares), (b) Drude weight (triangles) 
and kinetic energy per site (circles) for the FKL model for 
L = 10 cluster, 6 conduction electrons, as a function of Jn/t. 
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FIG. 9: (Color online) Optical conductivity (in units of e^) 
for the 2 X 10 cluster with 4 holes, ei = —1, (a) t = 0.4, (b) 
t = 0.9, and for several values of J indicated in the plot, (c) 
Optical conductivity for the ID extended t — J model on the 
10-site chain, J = l, ei = — 1,4 holes, and several values 
oi t/ J as indicated in the plot, (d) Optical conductivity for 
the FKL model (Eq. [S} on the 10-site chain, 6 conduction 
electrons, and several values of Jn/t indicated on the plot. 



maximum starting at At. For a fixed this maximum 
first increases as J is increased from 0.2 to 0.4, and then 
it decreases as J is further increased. The position of this 
maximum is shifted to higher values of lj. In Fig. [DJc), 
cr(w) is shown for the ID i — J model, on the 10-site chain 
with 4 holes. Notice that J is actually K in the notation 
of the Hamiltonian ([1]) and hence we adopted the value 
J = 1. Also, for the sake of comparison, an on-site poten- 
tial every two sites, with ei = —1, was added to the usual 
t — J Hamiltonian. The behavior of cr(a;) in both models 
is clearly different. In the extended t — J model there 
is a peak located at w ^ 2i instead of the broad feature 



observed for the FTJ model. In the present model, u{uj) 
is more resemblant to the behavior in two-dimensional 
clusters. In two dimensions, the competition between 
kinetic and magnetic energies is more important than in 
one dimension, where in general spin-charge separation 
holds. 17(0;) for the FTJ model is also more similar to 
the one for the FKL model (O, as shown in Fig. [Hd) 
for an equivalent chain size and filling: the incoherent 
structure is broad, and the spectral weight below this 
region (but at finite lj) may correspond in both models 
to the presence of a pseudogap. However, in the FKL 
model there is a pronounced change of (t(w) when enter- 
ing in the FM region {Jh /t > 19) implying an important 
transfer of weight from the incoherent part to the Drude 
peak. This behavior in the FM region is consistent with 
that reported in previous studies ,—122,1^ and it has been 
understood within the DE models The small remaining 
incoherent part does not scale with J/z/t and it is almost 
certainly a finite size effect. This suppression of the inco- 
herent part of <j(w) is virtually absent in the frustrated 
t-J model for VOx. Notice also the much smaller scale 
of (j{u!) for the FTJ model with respect to those for the 
extended t — J and FKL models. These different scales 
are in part due to the fact that t was used as the scale 
of energy in the latter models while K was used in our 
model. However, even after correcting for this factor, 
there is still an order of magnitude between a^'^^{u;) in 
the FTJ model and the other two models. 

Finally, the Drude weight for the 2 x 10 cluster with 4 
holes, ei — —1, is shown in Fig. fTUf a) as a function of J, 
for t = 0.4 and t — 0.9. It is instructive to examine first 
the behavior of the integral / of a^'^^{uj) (second term 
in Eq. ^) which measures the weight of the incoherent 
part of the optical conductivity. This quantity, for a given 
value of t, presents a maximum just at the crossover be- 
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tween the S = and S > phases, and then it smoothly 
decreases upon further increase of J. However, this inte- 
gral is much smaller than the kinetic energy (notice the 
factor of 4 in Fig. [TO]). As a result, as it can be clearly seen 
in this Figure, the Drude weight qualitatively follows the 
behavior of the kinetic energy, shown in Fig. [Hb). That 
is, for a given value of t it presents sudden drops at the 
values of J at which the total spin of the ground state 
increases. The behavior of the Drude weight of model 
([1]) is quite different from that of the FKL model. In 
Fig. [TOT b) we show the Drude weight for the later model 
on the L = 10 chain with PBC, 6 conduction electrons, 
obtained by ED. Although the kinetic energy monotoni- 
cally decreases when Jn/t increases due to the competi- 
tion with the magnetic energy the Drude weight first 
decreases while still inside the IC phase and then it is 
strongly enhanced upon entering in the FM phase, ap- 
proximately at Jn/t ~ 19, as mentioned above. This 
sudden increase of the Drude weight is consistent with 
the noticeable suppression of the incoherent structure of 
the optical conductivity in the FM region, observed in 
Fig.Hd). 

V. CONCLUSIONS 

We have proposed a ID model to describe a mechanism 
suggested to explain the appearance of a FM phase in 
vanadate nanotubes upon doping. This model consists of 
a t- J chain coupled to a Heisenberg chain with frustrated 
exchange interactions. 

Magnetic properties show the presence of an ascending 
"staircase" of the total spin 5* as the magnetic couplings 
Ji — J2 = J increase for a fixed value of the hopping 
integral t. These phases with S > have a ferrimag- 
nctic origin and spin-spin correlations do not present the 
typical behavior of a ferromagnetic state but they cor- 
respond to an incommensurate phase which is described 
by a peak of the structure magnetic factor located at 
q = {{L±1)tt/ L, 0) independently of the doping fraction. 
Our results suggest that this fcrrimagnetic phase will sur- 
vive in the thermodynamic limit. We have also observed 



this staircase in S in the FKL model, a result which was 
not reported to our knowledge in previous studies on this 
modelJ^ii^ii^ However, we cannot exclude at this point 
the possibility of this result being a mere finite size effect. 
In this model, we have also observed that the magnetic 
peak moves from its IC value to its FM value in a step- 
like fashion as Jh / 1 increases at a fixed density. 

The optical conductivity for the ID model for VO^^ 
shows a broad incoherent structure with indications of 
the presence of a pseudogap, and spectral weight that is 
transferred to higher frequencies as J is increased for a 
given value of the hopping constant t. This behavior is 
different from the one observed in isolated t — J chains, 
and more relevant to the topic of ferromagnetic metals, 
it is also very different from the one observed in the ID 
FKL model. This different behavior in (7{uj) is translated 
to the Drude weight. In the FKL model there is a strong 
enhancement of the Drude weight upon entering in the 
FM phase, while in the present model the Drude weight 
essentially follows the behavior of the kinetic energy. 

In summary, in the proposed model for VO^;, at least 
in its ID version, there is not a true ferromagnetic or- 
der but rather an IC fcrrimagnetic one, and contrary to 
what happens in the FKL model, electronic transport is 
somewhat suppressed by this ferrimagnetic order. More 
realistic variants of the model here studied would include 
a three-chain ladder or a two-dimensional lattice, a Hub- 
bard interaction instead of the t-J one, and a dimeriza- 
tion on the t-J chain. By increasing the dimensionality it 
could be possible to achieve a ratio between Vi, V2 and 
V3 ions closer to the real system. It would be closer to the 
actual materials to include a direct hopping between V2 
sites. Even for the present model, we have only explored 
a small region in the parameter space, and in particular 
it seems promising to examine the case Ji ^ J2, as we 
plan to perform in future works. 
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